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The complexity of a quantum state may be closely related to the usefulness of the state for 
quantum computation. We discuss this link using the tree size of a multiqubit state, a complexity 
measure that has two noticeable (and, so far, unique) features: it is in principle computable, and 
non-trivial lower bounds can be obtained, hence identifying truly complex states. In this paper, we 
first review the definition of tree size, together with known results on the most complex three and 
four qubits states. Moving to the multiqubit case, we revisit a mathematical theorem for proving 
a lower bound on tree size that scales superpolynomially in the number of qubits. Next, states 
with superpolynomial tree size, the Immanant states, the Deutsch-Jozsa states, the Shor’s states 
and the subgroup states, are described. We show that the universal resource state for measurement 
based quantum computation, the 2D-cluster state, has superpolynomial tree size. Moreover, we 
show how the complexity of subgroup states and the 2D cluster state can be verified efficiently. 

The question of how tree size is related to the speed up achieved in quantum computation is also 
addressed. We show that superpolynomial tree size of the resource state is essential for measurement 
based quantum computation. The necessary role of large tree size in the circuit model of quantum 
computation is still a conjecture; and we prove a weaker version of the conjecture. 


I. INTRODUCTION 

While we all have a feeling of what “complex” means, 
it is notoriously hard to find quantitative measures. Fur¬ 
thermore, there are various types of complexity. The 
three main examples in present-day science seem to be 
computational complexity, process complexity, and state 
complexity. Computational complexity refers to the 
amount of resources required to perform a certain compu¬ 
tation task, be it in terms of time, memory space or num¬ 
ber of queries, contributing to different complexity classes 
in computer science [1]. Process complexity is often as¬ 
sociated with the chaotic (but not random) behavior of 
the process, interconnectivity of many components in the 
process, and possibly the phenomenon of emergence [2]. 
Finally, state complexity, the focus of this paper, refers 
to the amount of information that is required to describe, 
generate or simulate a state of a physical system. 

Why do we study the complexity of quantum states? 
From a foundational point of view, complexity could be a 
parameter to test the limits of quantum mechanics. The 
direct extension of quantum effects (coherent superposi¬ 
tions) to daily objects might result in bizarre paradoxes, 
as Schrodinger famously noticed. It is a current experi¬ 
mental trend to push the tests of quantum mechanics to¬ 
wards the macroscopic domain, see for example Refs. El- 
E]. In all these experiments, the superposition indeed 
involves large number of particles or excitations, never¬ 
theless the states produced are somewhat "simple”: some 
involve the superposition of a single degree of freedom, 
the center of mass; others target the GHZ state, or the 
Dicke state with few excitations, as ideal macroscopic 
states. The macroscopic objects of our daily experience 
are not only large in size, mass and number of particles, 
but at the same time also interconnected in a non-trivial 
manner: a cat, besides being large, is a complex object. 
Do complex objects still obey quantum physics? If yes. 


as most physicists would argue, can we create them in 
a controlled way? These questions loomed behind the 
discussion on the possibility of large-scale quantum com¬ 
putation. In order to refine this discussion, Aaronson [7] 
took a technical step and proposed the concept of tree 
size (TS), as a measure of complexity for quantum states. 
This highlights that, besides exploring the limits of quan¬ 
tum mechanics, quantum state complexity is a way of 
capturing the deep relation between complexity and com¬ 
putation. 

The origin of quantum speed up might be sought in 
some features of entanglement. Promisingly, early stud¬ 
ies showed that states used in various quantum algorithm 
display multipartite entanglement [5]; while states with 
little entanglement could be efficiently simulated with 
classical computing 13 HD]. Nonetheless, large entangle¬ 
ment is neither necessary m nor sufficient condition for 
quantum speed up: to the contrary, having too much 
entanglement might be detrimental to performing com¬ 
putation [12] ■ Measures of entanglement developed with 
other operational meanings do not seem to capture the 
computational power of the state. 

Another candidate is the phenomena of interference. 
Previous works |13j propose how to quantify interference 
with “ibits” and investigate how many ibits were “actu¬ 
ally used” in various quantum algorithms. The different 
amount of actually used ibits seems to explain the differ¬ 
ent amount of speed up in Shor’s and Grover’s algorithm. 
The relation with success probability in algorithm with 
imperfections was studied in M- 

In this work, we look into the relation between com¬ 
plexity of quantum states and quantum speed up. It 
seems intuitive that, in order to he useful in computa¬ 
tional tasks, a state must be complex to describe and yet 
be simple to prepare. Indeed, if on the one hand a state is 
simple to describe, it should be possible to simulate it ef¬ 
ficiently with classical computers; on the other hand, the 
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preparation of the state from easily available resources is 
part of the overall computation process. Here, we focus 
on the first aspect: how to quantify the complexity of 
describing a quantum state! 

Among the different measures of state complexity, 
quantum Kolmogorov complexity is defined by length of 
the shortest possible program that would generate the 
state [IMIS]- This very common definition suffers from 
the setback that it is not computable. Moreover, Kol¬ 
mogorov complexity captures the complexity of gener¬ 
ating the state. The tree size (TS) complexity that we 
mentioned before and that we are about to discuss re¬ 
lates more to the description and simulation complexity 
of quantum states. The most common way to represent 
quantum states is the Dirac notation. Tree size complex¬ 
ity can be understood as the size of the minimal descrip¬ 
tion using this notation. 

This article provides a concise summary of our knowl¬ 
edge of the tree size, as well as some new results on com¬ 
plex states with superpolynomial tree size, verification 
of complex states, and the connection between tree size 
complexity and the power of quantum computation. The 
definition of tree size complexity is given in Sec. 2. Next, 
we review the works on states of two, three and four 
qubits, and describe the most complex states according 
to this measure. Moving to the case of n qubits, we first 
discuss a few examples of simple states with polynomial 
tree size. In Sec. 4, a theorem by Raz for showing su¬ 
perpolynomial lower bound on multilinear formula sizes, 
which in turn lower bounds tree size, is revisited. With 
this theorem at hand, we show some families of states 
with superpolynomial tree size: the Immanant states, the 
Deutsch-Jozsa states and the subgroup states. Based on 
numerical evidence, we construct an explicit example of a 
subgroup state with superpolynomial tree size. More im¬ 
portantly, the tree size of the 2D cluster state is shown 
to be superpolynomial. In Sec. 5, we describe how to 
verify the superpolynomial tree size of the complex sub¬ 
group states and the 2D cluster state with polynomial 
effort via measuring a witness. The possible relation be¬ 
tween state complexity and quantum computation speed 
up is discussed in Sec. 6. Finally, we offer a list of open 
problems and technical conjectures. 

II. TREE SIZE 

A. Definition and basic properties 

Just as bits to classical information, qubits are the ba¬ 
sic building blocks of quantum information. Any n-qubit 
pure states can be written in the computational basis 
with at most 2” coefficients: 

where each n x 1 vector x in {0,1}" is identified as a bit 
string in |a;). This decomposition on a computational ba- 
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Figure 1. A rooted tree of a three-qubit biseparable state. 


sis is not the most compact way of writing a pure state. 
In the case of two qubits, the most economic represen¬ 
tation is given by the Schmidt decomposition. This de¬ 
composition can be iterated to deal with multi-partite 
states, but already for three qubits a different ad hoc de¬ 
composition is more compact (Ej. The shortest possible 
representation of a multiqubit state in Dirac notation is 
given by the minimal tree size introduced by Aaronson 
in Ref. [7| as a measure of complexity of a pure state: 
Any multiqubit state written in Dirac’s notation can be 
described by a rooted tree of 0 and -I- gates; each leaf 
vertex is labelled with a single qubit state a |0) -I- ,5 |1), 
which needs not be normalized. The three-qubit bisep¬ 
arable state |0) (|00) -I- |11)), for instance, is represented 
by the tree in Fig. The size of a tree is defined as the 
number of its leaves: thus, the size of the tree of Fig. 
is five. A given state can have many different tree rep¬ 
resentations (for instance, the biseparable state in Fig. 
can also be written as |000) -I- |011) whose size would be 
six). The minimal tree of a state is the tree with the 
smallest size that describes it; and the tree size TS{|V')) 
is the size of the minimal tree. 

This measure of complexity is in principle computable, 
though we lack efficient algorithms. Moreover, a rela¬ 
tion with multilinear formulas leads to lower bounds on 
the tree size. It is thus possible to show that the tree 
size of some states is definitely superpolynomial in the 
number of qubits n |7| [50] . Such states can be consid¬ 
ered genuinely complex, in the sense that they cannot 
have a polynomial (i.e. computationally efficient) repre¬ 
sentation in Dirac notation — nor in a matrix-product 
representation with matrices of constant size . In con¬ 
trast, if the TS scales polynomially with n, the state can 
be considered simple as it can be described efficiently on 
a classical computer. Before moving to the n-qubit case, 
let us familiarize ourselves with TS by looking at the 
states of a few qubits. 
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B. Tree size of few-qubit states 

One observation that is very useful for finding TS of 
few-qubit states is the fact that TS is invariant under 
invertible local operations (ILOs). Formally [20] . 

Proposition!. If \ifj) = Ai ^ ^ An {(j)) , where all the 

single-quhit operators Ai are invertible, then TS(|'!/')) = 
TS(|<^)). 

Any two states that can be transformed to each 
other by ILOs as above are said to be equivalent under 
stochastic local operation and classical communication 
(SLOCC). The above proposition implies that all states 
in a SLOCC equivalent class have the same TS. 

a. Two qubits: Any two-qubit state can be written 
in the Schmidt decomposition as m- 

|^) = c| 0 )®| 0 ) + s|l)®|l), ( 2 ) 

where c and s are nonnegative real numbers satisfying 

-I- = 1, and {|0) , |1)} form an orthonormal basis. 

The state is said to be separable if one of the coefficients 
c or s vanishes and entangled otherwise. The Schmidt 
decomposition has size at most 4, hence the TS of any 
two-qubit state is at most 4. There are only two different 
rooted trees of size at most 4 that can describe a two- 
qubit state, which are shown in Fig. From this figure 
we see that a two-qubit state has TS = 4 if it is entangled 
and TS = 2 if it is separable. This concludes the case for 
two qubits. 

b. Three qubits: For three qubits, a useful decompo¬ 
sition that has a similar role as the Schmidt decomposi¬ 
tion does for two qubits is the canonical form derived by 
Acin et al. [E]: Any three-qubit state can be written as 

|'0) = COS0 |000) -I- sin0 |1) (cosa; |0'0") -I- sinw |l'l")|5) 

where the prime and double prime indicate different 
bases. The TS is upper bounded by the size of this de¬ 
composition and thus is at most 8 . Similarly to the case 
of two-qubit states, we first find all the possible trees for 
three qubits with at most eight leaves, and then try to 
see which one is the minimal tree of a given state. 

As stated in Proposition since all the states in a 
SLOCC class has the same TS, we need only find TS of 


Figure 3. Minimal tree of the (unormalized) |W) state 

one state in a class to know TS of all the states in that 
class. For three qubits, it is known that the pure states 
can be categorized into six different classes: the product 
class, three biseparable classes due to permutation, the 
GHZ class and the W class [32]. Examples of states in 
these classes are, respectively, 

|P) = 1000), 

|B)=^|0)(|01) + |10)), 

|GHZ) =^|000) + |111), 

|W) (1001) + 1010)-fIlOO)). (4) 

So, a state \if) is said to be in a particular SLOCC class, 
say the W class, if there exist ILOs Ai, A 2 , A 3 such that 
\iIj) = Al O A 2 ® A 3 |W). 

An exhaustive search [35] shows that TS(|W)) = 8 , 
TS(|GHZ)) = 6 , TS(|B)) = 5, and TS(|P)) = 3. So 
the TS of a three-qubit state can adopt only one of these 
four different values depending on what SLOCC class the 
state belongs to. 

The most complex three-qubit states are obviously the 
ones in the W class, whose minimal tree is drawn in 
Fig. [31 Interestingly, this complexity class is unstable 
in the sense that an arbitrarily small deviation could 
bring the W state to a state in the GHZ class [33| [24] . 
For studying how TS changes in the presence of fluctu¬ 
ation, we define a smoothed version of the tree size over 
a small neighbourhood of the desired state. For a posi¬ 
tive constant e < 1 , the e-approximate tree size of \ip) is 
the minimal tree size over all pure states \4i) such that 

I (</^’|^) P > 1 — e, that is, 

TS,(|^))= min TSm). (5) 

l( 0 lb>r>i-e 

Since in practice we cannot know a state with arbitrary 
precision, the TS^ is a more physical measure. The in¬ 
stability of the |W) state could be now phrased precisely 
as follows: For e arbitrarily small, there exists a state 
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I'/'GHz) in the GHZ class such that | (<^ghz|W) P > 1 — e. 
Therefore, TSe(|W)) = 6 . 

c. Four qubits: As in the case of three qubits, 
SLOCC equivalent classes can be used to find the tree 
size of four-qubit states. In Ref. [53] it is shown that the 
maximal TS of four-qubit state is 16 and a set of criteria 
for identifying whether a given state has this maximal 
TS is given. The class of most complex four-qubit states 
were found to belong to a SLOCC class not described in 
previous inductive classihcations [5^. These states have 
the following minimal decomposition: 

1^) = kl2) \Fm) + W2a) . (6) 

where \(f)) and \q)) are two-qubit entangled states. Noting 
the switching of order of qubits in the second branch, this 
form seems to preclude a recursive construction of the 
most economic description in terms of tree size. Forms 
that look recursive do require 18 leaves for some states. 

An example of the most complex four-qubit states with 
TS = 16 is 


S C B C W C GHZ. From the definition of TS for mixed 
states one sees that states belongs to the class S have 
tree size 3, B \ S tree size 5. For states in GHZ \ B, the 
tree size is 8 if a W state is required in the decomposi¬ 
tion, otherwise it is 6 . As an example we look at a family 
of one-parameter mixed state, the so called generalized 
Werner state m, 

p(p)=p|GHZ)(GHZ| + (l-p)^. (9) 

By looking at what SLOCC class p{p) belongs to for dif¬ 
ferent values of p, it is shown that TS = 3 (p G S) when 
p < 1/5, TS = 5 (p G S \ B) when 1/5 < p < 3/7. With 
an obvious decomposition into GHZ and product states, 
TS = 6 for 3/7 < p < 1, even though p G W \ B when 
3/7 < p < Pw and p G GHZ \ W when p > p^,, where 
p„ « 0.695 5427 [23]. 

III. SIMPLE STATES 





(10110) -b 10101) -b |1001) -b |1010)) 
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In this computational basis expansion its size is 24, but 
it can be shown that its minimal decomposition has in¬ 
deed the form given in Eq. ([^ with size 16. This state 
has already been created in experiments with four-photon 
down conversion [ST] [28] . 

When fluctuation is taken into account, the maxi¬ 
mal e-tree size of four-qubit states reduces to 14 for 
0 < e < 1/12: Any state in the most complex class can 
be e-approximated by a state with has the decomposi¬ 
tion |0) |GHZ) -b |1) |GHZ') where GHZ and GHZ' are 
two states in the GHZ class of three-qubit states. 

d. Mixed states: The concept of tree size may be 
extended to a mixed state p = IV'j) (V'jl as follows: 

TS(p) = min{maxTS('(/i)}, ( 8 ) 

i 

where the minimization is done over all the possible pure 
state decomposition of the mixed state p. This is the 
same approach for extending an entanglement measure to 
mixed states discussed in Ref. [29]. The intuition behind 
this definition is that the tree size of a mixed state should 
be at least as complex as the most complex pure state in 
its decomposition. 

According to the classihcation of three-qubit mixed 
sates introduced by Acm et al. [ 30 ], there are four differ¬ 
ent SLOCC classes: Class S, the set of states that can be 
written as combination of pure separable states; Class B, 
for states that can be written as combination of separa¬ 
ble and biseparable states; Class W, for states that can 
be written as combination of separable, biseparable and 
W states; and class GHZ, for states that can be written 
as combination of all possible three qubit states. Clearly, 


We now move to the case of n-qubit states and con¬ 
sider how TS scales as n —>■ oo. A state (more precisely, 
a family of states indexed by n) is simple if its TS scales 
polynomially with the number of qubits n. For showing 
that a state is simple, it suffices to find an explicit decom¬ 
position with polynomial size. Some examples of simple 
states are given in Table [l] 


Table I. Summary of n-qubit simple states 


Product state 
GHZ„ states 
Dicke states 
MPS with rank y 


n 

2 n 


0{n^) 

0(n‘°s2x) 


The product state |0)®" is the simplest state in terms 
of tree size. It is usually regarded as the input for 
the circuit model of quantum computation. Obviously 
TS(|0)®") = n, which is the minimal TS for n-qubit 
states. 

The n-qubit GHZ state, (|0)®” -I- |l)®”)/-\/2, which 
saturates most of the macroscopicity measures [35], has 
TS(|GHZ„)) = 2n, which is linear in the number of 
qubits. This is a clear evidence that complexity is a dif¬ 
ferent notion from macroscopicity. A maximally macro¬ 
scopic state can yet be very simple. 

The Dicke states \Dn^k) represents the equal superpo¬ 
sition of n-qubit string with k excitations; formally it 
is the (unnormailized) uniform superposition of all the 
n-bit strings with Hamming weight k: 

|i5„,,) = ^(g)|0),^{„j(g)|l),,{„j (10) 

fa} 

where the summation is over {a}, all the distinct k ele¬ 
ment subset of { 1 , • • • , n}. 
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We show that |-Dn,fc) has tree size O(n^). To see this, 
one can consider the uniform superposition of the follow¬ 
ing Fourier form (omitting normalization): 


fc-i 

\ilJn,k) = X! ( I®) +exp( 

j=o 


2TTij , 


ll) 


0n 


( 11 ) 


with tree size = 0{kn). A direct expansion 

yields {\Dn,p) Ei=d exp(27rijf . When p = mk, 

for some integer m, exp(27rij|) = 1; when p is not a 
multiple of k, exp(27rzj|) = 0. Hence, \ilJn,k) = 

J2m=’o \^n,mk)- For k > n/2, m can be only 0 or 1 , thus 
\Dn,k) = |'0n,fe) ~ |0)®”. For k = m can be 0, 1 and 
2; thus \Dn,k) = \^n,k) - |0)®” - For k < by 

interchanging 0 and 1 we obtain the fc > f case. So for 
any k, TS(|D„_fc)) = 0{v?). The n-qubit W state, which 
is \Dnp), though representing the most complex class for 
the three-qubit case, has polynomial tree size 0{n?). 

Finally, Matrix Product States (MPS) are a well- 
studied family because they approximate well the ground 
state of one-dimensional gapped Hamiltonians [5H1IM] . 
The tree size of an MPS is related to the bond dimen¬ 
sion X- -A. recursive argument provides the upper bound 
for the TS of an MPS [20]: consider the fol¬ 
lowing form of an MPS: 


IV')=(g)(4^|0)+A«|l)), (12) 

where and are matrices of dimension at most y. 
By partitioning the qubits into two halves, we have 

IC/ 2 ) IC/ 2 ) > (13) 

S=1 

where now MPS of n/2 qubits. We can see 

that TSd'i/n)) < 2x ■ TS(|'(/)„/ 2 )). By repeating this par¬ 
titioning, we have TSdV'n)) = 0 (( 2 x)*°®") = 0 (n'°s^^). 
Thus, if X is bounded as n increases, then TS is poly¬ 
nomial. One example of MPS, the ID cluster state, has 
X = 2, hence its tree size is O(n^). Note that the same 
recursive argument applied to a more general form of ten¬ 
sor network states, the projected entangled pair states 
(PEPS), gives the superpolynomial upper bound 
and for 2D and 3D PEPS respectively; and in¬ 

deed, as we are going to discuss later in Sec. |VlJ PEPS 
that are universal for measurement-based quantum com¬ 
putation (MBQC) should have superpolynomial tree size, 
assuming that factoring is not in P. 


IV. COMPLEX STATES 
A. Methods to obtain lower bounds on tree size 


One of the advantages of tree size as a complexity mea¬ 
sure is that there are tools for proving lower bound on 
tree size, hence certifying complex states. One way is 
to use counting argument as Aaronson did in Theorem 7 
of [7] . The fact that there are fewer states that has poly¬ 
nomial tree size than there are in the whole Hilbert space 
(or the state space of interest), some states are bound to 
have superpolynomial or even exponential tree size. 

Another method, which will be discussed more often 
in this paper, is a theorem first proved by Raz in the 
context of multilinear formula size (MFS) [71135]. Al¬ 
though counting arguments could show that states with 
superpolynomial tree size must exist, but Raz’s theorem 
allows us to construct explicit examples. Let us present 
here this important theorem, first in Raz’s original for¬ 
mulation, then in an equivalent way in terms of Schmidt 
rank. 

A multilinear formula is a formula that is linear in 
all of its inputs. The MFS of a multilinear formula is 
defined as the number of leaves in its minimal tree rep¬ 
resentation similar to the tree size of a quantum state. 
Consider a multilinear formula / : {0,1}" —>■ C, let P 
be a bipartition of the input variables cci,-- - ,Xn into 
two sets, j/i,--- , 2/„/2 and Zi,--- tZ^ii- We now view 
f{x) as a function fp{y,z) : {0,1}"/^ x {0,1}"/^ —>■ C. 
Then denote by Mf^p the 2"/^ x 2”/^ matrix whose rows 
and columns are labeled by y and 0 S { 0 , 1 }"/^, re¬ 
spectively. The entry (y, z) of this matrix is defined as 
Mf\p{y,z) = fp{y,z). Finally, let rank(M^|p) be the 
rank of Mf^p over the complex numbers, and P be the 
uniform distribution over all the possible bipartitions P. 
We now state Raz’s theorem [55] : 


Theorem 2. If 


Pr 

PGP 


rank(M^|p) >2 2 




(14) 


then MFS(/) = 

For any quantum state \tjj), we can define the associ¬ 
ated multilinear formula /^(a:) = {x\ip). Note that this 
formula computes the coefficients in the computational 
basis expansion of Given a tree representation of a 
quantum state, a tree for the associated multilinear for¬ 
mula can be obtained by interchanging | 0 )j —>■ {1 — Xi) 
and |l)j —>■ Xi- Thus, given the minimal tree of quantum 
state, we can obtain a tree for the associated multilinear 
formula with the same size. The true MFS of the formula 
can only be smaller, therefore [ZlEO]: 

Theorem 3. MFS(/^) < TSdz/)). Therefore, if fi/, sat¬ 
isfies Raz’s theorem, then TS(|^)) = 

In fact, in the original paper of Aaronson [7], he showed 
that the other way of the inequality is also true up to n, 
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TS(|^/>)) = 0(MFS(/^) + n). But for the purpose of the 
rest of this article, MFS(/^) < TS(|^/>)) is sufficient. 

For its application in complexity of quantum states, we 
rephrase Raz’s theorem in terms of the Schmidt rank, a 
well-known concept in quantum information | 21 j : 


Theorem 4. For a pure quantum state of n qubits, con¬ 
sider all the uniformly distributed (§, §) bipartitions, if 


Pr 


T.-r.l/S 

SR > 2^^ 




(15) 


where SR is the Schmidt rank of a particular partition, 
then TS = 


The statement follows indeed from Raz’s theorem, be¬ 
cause partitioning of the input x of the associated multi¬ 
linear formula is the same as partitioning the qubits 
of the state l"^). Note that Mf^\p is a matrix each of 
whose element is a coefficient of the state IV') in its com¬ 
putational basis. Thus, the rank of Mf^\p is exactly the 
Schmidt rank of the state {ip) for the bipartition P |21) . 

Interestingly, from the point of view of complex sys¬ 
tems and statistical physics, multipartite entanglement 
were found to be related to the average entanglement 
across equal bipartitions [? ? ? ? ]. Instead of 

the Schmidt rank, the distribution of purity of the par¬ 
tial states over all equal bipartition was studied in those 
works. This suggests a possible deeper link between tree 
size complexity and multipartite entanglement. 

With the help of these theorems, we shall identify some 
explicit multiqubit states with superpolynomial tree size. 


B. Immanant states 


An explicit family of states with superpolynomial tree 
size can be constructed based on the immanant of a 
(0,1) matrix [^. Consider the case when the number 
of qubit is a square number, n = mf, for each bit string 
|a;) = |a;i,...,x„) we arrange the bits xi,... ,Xn row by 
row to an TO X m matrix M{x) such that 


M(x) 


/ Xi 
^m +1 


X2 

Xm+2 


Xm 

X2m 


(16) 


where a is an element of the symmetric group Sm of 
all the to! permutation of { 1 , 2 ,...,to}, and c^ is the 
corresponding complex coefficient. When Co- = 1 for all a 
the immanant reduces to the permanent, and when = 
1 for even permutations and —1 for odd permutations it 
reduces to the determinant. It is proved in Ref. [20] that 

Theorem 5. The Immanant states as defined above have 
TS = if the coefficients c,, are all nonzero. 

The proof relies on Raz’s technique [35] to show 
that the multilinear formula size of the immanant with 
nonzero coefficients is superpolynomial, and this theorem 
follows immediately from Theorem]^ 

The permanent and the determinant states, 

|Perm„) = ExGfo,!}’* Perm(M(x)) \x) , 

|Det„) =ExG{o,i}"Det(^(2;))|2^)> (19) 

are two examples in this family of complex states. 

The smallest known formula for computing permanent 
is the Ryser’s formula [? ], which is multilinear: let S be 
one of the 2 ’" subsets of ( 1 , 2 ,..., to} and [S'] the number 
of its elements, then the permanent of the matrix M is 

m 

Perm(M) = E HE My. ( 20 ) 

s i=ijes 

By substituting this formula to the permanent state and 
carrying out the summation over x we obtain a decom¬ 
position with size n^2^. We conjecture that the tree 
size of the permanent state is TS(|Perm„)) = see 

Sec. |VIIB0g| for detailed discussion. 

A common confusion sometimes arises: why do we 
treat the permanent state and the determinant state on 
the same footing while determinant is known to be much 
easier to compute than the permanent: In fact, there ex¬ 
ists a formula that computes the determinant of a to x to 
matrix with size O(to^) [37]. This does not contradict 
with Raz’s result of since the optimal algo¬ 

rithm does not use a multilinear formula; and only mul¬ 
tilinear formulas can be used to find an upper bound on 
the TS of the corresponding state. 


C. Deutsch-Jozsa states 


\fXn—m Xn—m-\-l ’ ‘ ’ Xn j 

so M(x)ij = Xm(i-i)+j- The immanant states are defined 
in its computational basis expansion as 

|Immn) = ^ Imm(M(x)) |x) . (17) 

xGfO,!}" 

Here the immanant Imm(M) of a matrix M is given by 

m 

Imm(M) = E (18) 

(7GSm i —1 


The Deutsch-Jozsa algorithm outperforms its classical 
counterparts in the deterministic case [2T]. It is an algo¬ 
rithm that solves the following hypothetical question: A 
function / : { 0 , 1 }" —)■ ( 0 , 1 } is called balanced if exactly 
half of its input is mapped to 0 and the other half to 1 , 
and constant if all the inputs are mapped to 0 or 1. Given 
the promise that the function is either balanced or con¬ 
stant, how many queries do we need to find out whether 
the function is balanced or constant? Classically, in the 
deterministic and worse case scenario, it requires 2 "“^-|-l 
queries, in which case the function outputs all 0 or all 1 
for the first 2 ”“^ queries. 
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The Deutsch-Jozsa algorithm solves the quantum ver¬ 
sion of this problem with only one query, which is ex¬ 
ponentially faster than the classical algorithm. In the 
quantum version, a query is replaced by the quantum 
oracle |a;) |?/) —>■ \x) \y (B f {x)) . In this algorithm, one 
first prepares the input state as lO)" | 1 ), then applies the 
Hadamard transformation to all the registers, resulting in 
2 -("+i )/2 |a;) (|0) - |I)). After applying the or¬ 
acle, the state becomes X]a:G{o i}" 1 ^) (l/(^)) ~ 

|I 0 f{x))). Since f{x) is either 0 or 1, we can simplify 
this to IV') = li’Dj) 0 |—), where 

\^oj) = ^ E (21) 

sGlO.l}" 

The last qubit register can be left out at this 
point. Applying the Hadamard transforma¬ 

tion to all the qubits once again, we have 

2 -"E,g{ 0 .i}" where x ■ y 

represents the sum of bitwise product. Finally, a projec¬ 
tion onto | 0 )” has probability | 2 “"X)a;e{o i}’*(“ 1 )'^^^^ 
which evaluates to 1 if f{x) is constant and 0 if f{x) is 
balanced. This concludes the algorithm, now we switch 
the focus to the tree size of the state \ipDj)- 

If / is constant, then the state {tpaj) is a simple prod¬ 
uct state 1+)" with tree size n. If / is balanced, we would 
like to show that an overwhelmingly large fraction of bal¬ 
anced functions correspond to states with superpolyno¬ 
mial tree size. Consider a function / randomly drawn 
from the uniform distribution of all the ( 2 ^/ 2 ) balanced 
functions, let P be a random equal bipartition of the in¬ 
put X into y and z, then Mj|p the 2”/^ x 2"/^ matrix 
whose entries are Note that for a balanced 

/, the matrix Mf^p has exactly half entries equal to 0 l 
and half equal to —1. Let £i be the event that Mj^p 
has full rank 2"/^, we need to compute the probability 
that £i happens, in order to see whether the balanced 
function / leads to a state with superpolynomial TS (c.f. 
Theorem . 

Let us call a matrix with exactly half entries equal to 
1 and the other half —I a balanced (1,-1) matrix. Denote 
by Mp a random balanced (1,-1) matrix, Mp can be 
chosen by first drawing a random balanced function /, 
then picking a random bipartition P and assigning Mp = 
Mf\p. Now let £2 be the event that Mp has full rank, 
we have 

PTi£2)=Y.PT{f)FTi£,\f). ( 22 ) 

/ 

Next, we split the set of / into those which give rise to 
a complex state (i.e. satisfy Raz’s theorem) and those 
which do not. Explicitly, let C = {/| Pr(fi |/) > g}, 
where g is a constant to be specified later, and C be 
the complement of C\ then 

Pt{£ 2 ) = Y. Pr(/) Pr(£i|/) 0 ^ Pr(/) Pr(£i |/). (23) 
C C 



Figure 4. The probability of the event £2 versus the number 
of input bits n. £2 is the event that a random 2"'^^ x 2 "A 
balanced (1,-1) matrix has full rank. 


Since Pr(£’i|/) < 1 for all / and Pr(£’i|/) < q for all 
f G C, we have 

Pr(52)<5]Pr(/)+g^Pr(/). (24) 

C C 


Note that the sums of the probability that / is cho¬ 
sen from C and C give the fraction of states in the 
respective sets, that is, P^(/) = ^ and 

^gPr(/) = ^ ~ ^C- Substituting this to 

the above inequality, we arrive at 


Fc > 


Pr(g2) - q 

1-g 


(25) 


Thus, to know how large Fq is we need to know Pi{£ 2 ), 
the probability that Mp is invertible where Mp is a ran¬ 
dom 2"/^ X 2”/^ balanced (1,-1) matrix. Our numeri¬ 
cal evidence shows that Pr(£ 2 ) approaches 1 quickly as 
n becomes large (see Figure [^. If one believes that 
Pr(f 2 ) ~ 1 for large n, which is strongly suggested by 
the numerical evidence, then by setting g to a constant 
not close to 1, say 0.5, we see that Fc « 1. This means 
that nearly all balanced functions give rise to states with 
superpolynomial tree size. 

One may argue that the large tree size that arises from 
the Deutsch-Jozsa algorithm has its root in the oracle’s 
access to completely-random balanced function. The link 
between large tree size and the usefulness of the algo¬ 
rithm is unclear. Nonetheless, this provides us with an 
example of complex states that appear in a quantum al¬ 
gorithm. More on the relation between state complexity 
and quantum computation will be discussed in Sec. lYll 


D. Shor’s states 

Shor’s algorithm factors an integer N in time 
0((logAf)^) [38], which is exponentially faster than the 
most efficient known classical algorithm [33]. Do states 
arising from this algorithm have superpolynomial tree 









size? Aaronson showed the answer is yes assuming a 
number-theoretic conjecture [7j. To factorize N, pick a 
pseudo random integer s < N, coprime to iV, consider 
the Shor’s state of n = log(A^) qubits, which is given by 


1 2--1 

^ I*”) 1^’’ 

r—0 

To simplify the proof of the lower bound on the TS of 
the Shor’s state, it is convenient to measure the second 
register. Since a measurement in the computational basis 
does not increase tree size for any outcome of the mea¬ 
surement, we can assume that the measurement outcome 
to be 1. Then, the state of the first register has the form 

\p^) = 


where p is the order of s modulo N and I = [( 2 " — 1) /p\. 
Here pi is represented in binary with n bits, so |pZ) is a 
n-qubit state. TS(|pZ)) provides a lower bound for the 
tree size of the state of the two registers given in ( [ 26 | . 

The associated formula for this state is a function of 
a n-bit string such that fn,p{x) = 1 if a: = 0 mod p 
and fn,p{x) = 0 otherwise. MFS(/„_p) lower bounds 
TS(|pZ)), so we shall focus on this formula. Aaron¬ 
son showed MFS(/„,p) = assuming the fol¬ 

lowing number-theoretic conjecture [7]: there exist con¬ 
stants 7,(5 G ( 0 , 1 ) and a prime p = H( 2 "’ ) for which 
the following holds. Let the set A consists of ele¬ 
ments of { 2 °, • • • , 2 "“^} chosen uniformly randomly. Let 

S consists of all 2 ”^ sums of subsets of A, and let S 
mod p = {x mod p : x G S}. Then 


Pr 

A 


[S' mod p| > (1 -I- 7 ) 


2J 


-o(logn) 


(28) 


E. Subgroup states 


then at least 4% of these matrices give rise to subgroup 
states with superpolynomial TS. 

Let us describe briefly how to prove that a subgroup 
state has superpolynomial tree size. Consider a ran¬ 
dom equal bipartition of a; = {xi,... ,Xn} into y = 
{yi, ■ ■ • , 2 / 71 / 2 } and 0 = {zi, ... ,z„/ 2 }. Denote by Ay the 
n/2 X n/2 submatrix of the columns in A that applies 
to y (see Eq. @), and similarly Az the submatrix of 
the columns that applies to z. Then, the element of the 
partial derivative matrix Mf\p(y, z) is 1 when Ax = 0 
mod 2 , which means Ayy + AzZ = 0 mod 2 , and 0 oth¬ 
erwise. So long as both Ay and A^ are invertible, for 
each y there is only one unique value of z that gives 
Mf\p{y,z) = 1. In other words, Mf\p is a permutation 
of the identity matrix, hence it has full rank. Based on 
this observation, one sees that 

Theorem 6 . Let A be a n 12'An binary matrix and S = 
ker(A) over the field Z 2 . For random equal bipartitions 
of X into y and z as described above, if Ay and Az are 
both invertible with probability then TS(|S')) = 

n^(i°g”). Moreover, TSedS”)) = with e < 1 —/i„, 

where pn = 2 “^"'/^^ 

Proof. The first part follows from the fact that both Ay 
and Az being invertible implies that A/^|p has full rank. 
If this happens with probability then Raz’s 

theorem is satisfied, hence TSdS*)) = 

For the second part, we use a lemma proved by Aaron¬ 
son in Ref. [7]: Denote by \ip) a state close to a complex 
state \S) that satisfies theorem]^ such that | (iplS) ^ < 
1 — e. Then, for a fraction of of all equal bipar¬ 

titions, the rank of the partial derivative matrix is 

rank(My,|p) > (1 — e)2"^^. (31) 

In order to satisfy Raz’s theorem, we want rank(M^|p) > 
2 n/ 2 -(n/ 2 )i/®/ 2 _ ^ comparison with the above equation 
gives e < 1 — pn where pn = 2“^”/^)^'^*/^. Therefore, 
TSflfj}) = if e < 1 - 2 -("/ 2 )'^V 2 , q 


Subgroup states used in quantum error correction also 
exhibit superpolynomial tree size. Let the element of Titf 
be labeled by n-bit strings. Given a subgroup S' C Z 2 , a 
subgroup state is defined as 

One way to construct a subgroup state is by consider¬ 
ing the subgroup to be the null space of a ( 0 , 1 ) matrix 
over the field Z 2 . Given a n/2 x n binary matrix A, a bit 
string X is in the null space of A if 

Ax = 0 mod 2 ; (30) 

and the subgroup state is the equal superposition of all 
such bit strings. Aaronson shows in Ref. [7] that, if A is 
drawn from the set of all possible n /2 x n binary matrices. 


Since pn is exponentially small in n^/®, one might think 
that most states in the Hilbert space satisfy | {ip\S) p > 
Pn, and hence Theoremcan be used to show that most 
states have superpolynomial TS. This is not correct: In¬ 
deed, if Itjj) is randomly and uniformly chosen from the 
Hilbert space according to a Haar measure, the probabil¬ 
ity that I (^/’|S') p > Pn is smaller than exp[—( 2 ” — l)/x„], 
which is exponentially small |40j . However, it is true 
that most states in the Hilbert space have exponential 
tree size, as showed by a counting argument in Ref. [7]. 

Aaronson hrst showed an explicit construction by Van¬ 
dermonde matrix that leads to a superpolynomial com¬ 
plex subgroup state [7]. Here we present a different con¬ 
struction of the matrix A, for which strong numerical evi¬ 
dence suggests that the corresponding subgroup state has 
superpolynomial TS. Gonsider the matrix Aj = (1|Q), 
where 1 is the identity matrix and Q a binary Jacob- 
sthal matrix, both of size n/2 x n/2. Jacobsthal matrices 
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Figure 5. The probability of both Ay and being invertible 
over random equal bipartitions of Aj. A,j is the qx2q matrix 
(1|Q), where Q is the Jacobsthal matrix of size q x q, where 
q = 3 mod 8 and is a prime. For large q, the probability 
approaches a constant around 30%. 


are used in the Paley construction of Hadamard matri¬ 
ces [41]. The binary version is defined as follows: For 
a prime number q, one can define the quadratic char¬ 
acter x(a) that indicates whether the finite field element 
a e Zq is a perfect square. We have x(a) = 1 if a = 6^ for 
some non-zero element & G Z^; and x(a) = 0 otherwise. 
Then j is equal to x(i — j). 

We study the partitioning of Aj into Ay and ran¬ 
domly. Numerical evidence (see Fig. shows that when 
q is a prime and q = 8k + 3 with fc S N, then , Ay and 
are both invertible with a probability approaching to 
a constant around 30%. From Theorem]^ we see that the 
subgroup state defined by Aj has TS = where 

n = 2q. 


F. 2D cluster state 
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Figure 6. A schematic diagram of measurement-based quan¬ 
tum computation. Starting from a 2D cluster state, single 
qubit measurements are performed. © represents a Z mea¬ 
surement, and other arrows represent measurements in the 
XY plane. The logical input state enters from the left and 
propagates to the right. Single qubit rotations and controlled 
gates are realized by a certain sequence of adaptive measure¬ 
ments, from left to right. Implementing a circuit on n qubits 
with m gates requires a cluster state of size 0(n)-by-0(m). 


stabilizer circuit is one that consists of only CNOTs, 7 r/2- 
phase gates and Hadamard gates. In the MBQC scheme, 
each of these gates can be implemented by measuring a 
constant number of qubits: 15 qubits for CNOT, and 5 
qubits for the phase gate and the Hadamard gate [43|. 
In order to obtain a n-qubit complex subgroup state, 
one needs to prepare a 0(n)-by-0(n^/logn) lattice (see 
Fig.[§, so the number of qubits in the 2D cluster state 
is = 0{n^ / \ogn). Since single qubit projective mea¬ 
surements only decrease tree size (c.f. the proof of theo¬ 
rem 


10 ), we have 


It is known that measurement-based quantum com¬ 
putation (MBQC) on the 2D cluster state is as strong 
as the circuit model of quantum computation |421143j . 
In this scheme of computation, after the initial resource 
state is prepared, one only performs single qubit projec¬ 
tive measurements and feedforward the outcomes. The 
power of the computation seems to lie in the initial re¬ 
source state. Therefore, an initial state that is universal 
for quantum computation, such as the 2D cluster state, 
should be highly complex. It is conjectured in Ref. |7] 
that the 2D cluster state has superpolynomial TS. By 
studying the generation a complex subgroup state via 
MBQC on the 2D cluster state, we can prove that this 
conjecture is true: 

Theorem 7. The 2D cluster state of N qubits has TS = 

jYn(logAr)_ 


Proof. Suppose we aim to produce an n -qubit complex 
subgroup state |S'c) (as described in Sec. IVEI that has 
tree size These states are known to be stabi¬ 

lizer states [21]. Aaronson and Gottesmannshowed that 
any n-qubit stabilizer state can be prepared using a sta¬ 
bilizer circuit with 0(n^/log n) number of gates [H] . A 


TS(2D cluster) > TS(|5c)) = n^(^°8n) ^ p^niiog n'^^2) 


So, the Wqubit 2D cluster state has the superpolynomial 
tree size. □ 


We can also show that the e-tree size of the 2D cluster 
state is also superpolynomial: 

Theorem 8. For e < 1/2, TSe(2D cluster) = . 

Proof. Assume that we have prepared a state close 
to the 2D cluster state, |2Dj), such that the fidelity 
F(|2D),|2D,)) = |(2D|2D,)| > Then we ap¬ 

ply the same measurement sequence to the erroneous 2D 
cluster state as if we would to the ideal 2D cluster for 
preparing a complex subgroup state |5'c). Consider the 
state after one of the single-qubit measurement in the or¬ 
thonormal basis {|? 7 ), I??"'')}; the single-qubit projectors 
are Pq = \t]) {r]\ and Pi = lip^) (??■*■ |. We now show that 
one of these outcomes will increase the fidelity between 
the two cases. If the measurement outcome is not ob¬ 
served, the resulting states on the ideal and e-deviated 
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2D cluster states are: 

\2D) ^ p = Po \2D) {2D \Pj + Pi |2P) (2P| P/ 

= PO b) (??l ® IV’o) (V’ol +Pl (??-^| O IV’l) ('0l| , 

(33) 

|2D,) ^ a = Po |2P,) (2P,| Pj + Pi |2P,) (2D,| P/ 

= Po \rf) ( 77 I (g) IV’o) (t^ol +pi |r?-^) (r?-^| 0 |V’'i) (V’ll, 

(34) 

where |^/’o,i) and \'4’o^i) are the states of the remaining 
qubits in the cluster; and po,i and pp 1 are the probability 
of the measurement outcomes. Clearly, the above map is 
completely positive and trace-preserving (CPTP). Thus, 
the fidelity of these two states should not decrease due 
to monotonicity of fidelity under CPTP maps [2T], 


we verify that the produced state is superpolynomially 
complex? Full state tomography requires exponentially 
many operations and is hence not practical. Nonethe¬ 
less, for complex states that are stabilizer states, there 
exists a complexity witness that can he measured with 
only a polynomial number of basic operations. This wit¬ 
ness can be used for verifying the superpolynomial TS of 
pure states. Proving and verifying superpolynomial TS 
of mixed states remains an open problem. 

The subgroups states described in Sec. |IV E| belong to 
the class of stabilizer states. A n-qubit stabilizer state 
[S') is uniquely defined by n mutually commutative stabi¬ 
lizing operators in the Pauli group, gi,... ,gn, satisfying 
the eigenvalue equation: 

g. \S) = \S). (39) 


F(p,a)>F(|2D),|2D,)) = vT^. (35) 


With a bit of algebra, we can express F{p,a) in terms 
of the fidelity of the post-selected states for the same 
outcome: 

F{p,(j) = Tr-\/pCVpb^ 

= \/poPol (V'olV’o) I + y/piP'i\ I- (36) 

Let us denote x = max(| (^olf/’o) !> I (V'llV'i) I) to be the 
larger overlap between the two, then 

F{p, a) <x ( VpoPo -b VpiPi) < x, (37) 


since yjpop'n + \/pip\ < (po +Po +Pi+Pi)/2 = 1. Com¬ 
bining Eq. (36) and Eq. (37), we have 

a: = max (I (^/^olV’o) 1.1 (V'ilV''i) I) > Vl - e. (38) 


Therefore, for at least one of the outcomes, we have a 
non-decreasing fidelity on the unmeasured parts of the 
states. For every measurements we post-select on the 
outcome that do not decrease the fidelity. Note that 
the complex subgroup states can be realized by a Clif¬ 
ford circuit, which can be implemented by a series of 
non-adaptive measurements. This means that, regard¬ 
less of the outcome, the state obtained from the ideal 
2D cluster is a complex subgroup state l^c) upto local 
Pauli operators. For the erroneous 2D cluster state, we 
would obtain a state IS'^) such that | {Se\Sc) \ > — £■ 

From theorem we see that when n is large enough, 
TSdS”,)) = and hence TS(|2D,)) = iV^(iogW), 

if e < 1/2. Thus, TS,(|2D)) = for e < 1/2. □ 


The generators of the subgroups states can be read off 
from the corresponding matrix A. Let R = rank(A); then 
there are R linearly independent rows ri (1 < i < i?) in 
A. For the first R generators, one simply replace 0 by / 
and 1 by Z for each of the first R linear independent row. 
For example, if row is (0, 0,1,0), we write gi = IIZI, 
where the position of the operators denotes the qubit on 
which they operate on. The remaining generators can be 
found from the n — R linearly independent vectors Ci that 
span the null space of A. One replaces 0 with I and 1 
with X for each vector, and the generator is the ordered 
product of these operators. 

Proposition 9. The operators gi defined above are the 
generators of the stabilizer of \S). 

Proof. Recall that [S') is the uniform superposition of \x) 
where a; is a vector in the null space of A. For the first 
R generators, we have gi \x) = (—1)’'"'® |x) = |a;), for all 
X G ker(A), hence gi [S') = [S') for z = 1, • • • , R. For the 
generators obtained from the n — R linearly independent 
vectors Ci in the null space of A, we have gi\x) = |a; 0 c^), 
where 0 is the bitwise addition modulo 2. Note that Ci 
is in the null space of A, so ker(A) 0 = ker(A), and 

hence gi\S) = [S'). This shows that the gts stabilize [S'). 

For the commutation relation, it is obvious that the 
first R generators commutes with each other and so do 
the n—R obtained from the null space. It remains to show 
that gi from row Vi commutes with gj from Cj. ri ■ cj = 0 
mod 2 implies the number of positions where the entries 
of both ri and Cj are 1 must be even. The single-qubit 
operators in giPj at these positions are ZX = —XZ\ and 
since there are an even number of these pairs we see that 


9i9i ~ 9j9i- 


□ 


V. WITNESSING COMPLEX STATES 

In this section we address the problem of verifying the 
large TS of complex states. Suppose one wants to cre¬ 
ate complex states such as the complex subgroup states 
and the 2D cluster state in the lab, in reality the pro¬ 
duced states are at some distance away from the tar¬ 
get states due to experimental imperfection. How do 


Now we show how to construct a complexity witness 
based on the complex subgroup states. Consider a state 
\Sc) that satisfies Theorem]^ For large n, any n-qubit 
state 1^) such that | ('0|5'c) F >1/2 must have TS = 
n^(^°g"). The superpolynomial TS of these states can be 
verified by measuring the witness 

W = h - \Sc) {Sc\. (40) 
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A negative value of {W) implies that the overlap of the 
produced state and \Sc) is larger than 1 / 2 , and hence the 
TS of the produced state is superpolynomial. However, 
W as such is not measurable in practice, under the natu¬ 
ral constraint that only local measurements are feasible. 
If one decomposes W into a sum of locally measurable 
operators, the number of such measurements increases 
exponentially with the number of qubits [551 ESI ES]- 
Nonetheless, when \Sc) is a stabilizer state, it is possible 
to construct a stabilizer witness W with the following 
properties: If (Vb') < 0 then (H^) < 0; and W can be 
decomposed into a sum of a linear number of operators 
in the Pauli group, which in turn can be measured by a 
polynomial number of basic operations m- The stabi¬ 
lizer witness is defined as: 

n 

W = {n-l)t-J2g,. (41) 


To show that {W) < 0 implies {W) < 0, one considers 
all the eigenvalue equations of the form (39) but with 
possible eigenvalues ±1. This defines the set of 2" com¬ 
mon eigenstates of the generators piS. Since all the gen¬ 
erators are Hermitian operators, the common eigenstates 
are mutually orthogonal and form a complete basis. One 
can verify that, in this basis, the operator W'—2W is a di¬ 
agonal matrix with non-negative diagonal entries. Thus, 
W — 2W is a positive semi-definite operator; so {W’) < 0 
implies {W) < 0. If in an experiment the expectation 
value of the stabilizer witness W is found to be nega¬ 
tive, then one can certify that the produced state indeed 
has TS = 


While the witness W detects all complex states with a 
fidelity (with respect to \Sc)) larger than 1 / 2 , W' detects 
a smaller set. It is necessary to know how close to \Sc) a 
state lip) needs to be for {p;\ W |^/>) to be negative. If the 
required fidelity is exponentially close to 1 then no state 
would be detected by W in practice. For this purpose, 
we first expand \ip) as 


\^P) = ci\Sc) + C 2 \S^), (42) 


where is a state orthogonal to [S'c) and |cip-|-|c 2 p = 
1. We have 


{■>p\W'\->p) = n-l-n\ci\^ - \c 2 \'^^{S^\gi\S^). (43) 

i=l 

Since 1 pi is a positive semi-definite matrix, 
(5'-*-|gi > —1. Therefore, 

{'ip\ W' \'ip) <n — 1 — n|cip -I- n|c 2 p = 2n — 1 — 2n|cip. 

(44) 

Thus, (V’l W \'ip) < 0 when the overlap | {ip\Sc) P = 
|ci p > 1 — l/(2n). So, the loss of fidelity must be smaller 
than l/( 2 n) for a state to be detected by W. 

One needs to measure all the n generators to estimate 
{W). With the help of an ancilla qubit, all the gen¬ 
erators, each with two possible outcomes, can be mea¬ 
sured by applying a circuit of size O(n^) followed by 



Figure 7. A circuit for measuring all the generators giS (only 
two are shown here). The controlled-gi gate can be decom¬ 
posed into at most n two-qubit controlled-Pauli gates, and 
there are n generators to be measured. Projective measure¬ 
ments of the ancilla qubits in the computational basis give 
the outcome of gi. 


a measurement on the ancilla qubits [H] (see Fig. [^. 
These measurements need to be repeated to obtain the 
desired accuracy. Suppose the produced state has a fi¬ 
delity I ('0|5'c) P = 1 — a/( 2 n) with a < 1 is a constant, 
we have (W') < — (1 — a). If the random error in each 
Pi is 6g then 6W' = nSg. Thus, to be confident that 
{W) < 0 one needs nSg < 1 — a, or Jt; < (1 — a)/n, 
which is achievable with a polynomial poly(n) number 
of repetitions. Therefore, a correct negative expectation 
value of W can be obtained with polynomial effort. 

There is a similar stabilizer witness for detecting com¬ 
plex states close to the 2D cluster states. Indeed, the 2D 
cluster state has TS 1/2 = and is also a stabilizer 

state. Thus, the witness for the 2D cluster state has the 
same form as W, with the piS replaced by the generators 
of the 2D cluster state. These generators are described 
in Ref. [45] . 


VI. RELATION TO QUANTUM 
COMPUTATION 

One of the main motivation of this study is to inves¬ 
tigate the relation between state complexity and quan¬ 
tum computation. To elaborate on this, we can divide 
all the quantum states into four categories according to 
their preparation complexity and state complexity (see 
Fig. §. The set of states with large preparation com¬ 
plexity but small state complexity is presumably empty 
because preparing simple states should not be too dif¬ 
ficult. The states with small state complexity are not 
useful for quantum computation because they are too 
simple and hence a classical computer can simulate them 
efficiently. The states with large preparation complexity 
are not useful either because quantum computation with 
these states requires too much resource in space and time. 
States that are useful for quantum computation should 
be the ones that have large state complexity yet small 
preparation complexity. If tree size is a good measure 
of state complexity, then we might ask: is superpoly¬ 
nomial tree size a necessary condition for the state to 
provide advantage in some computational task? In this 
section, we are going to discuss this link in the frame- 






















12 



Figure 8. Dividing all quantum states into four categories ac¬ 
cording to their state complexity and preparation complexity. 
One out of the categories is presumably empty, two are not 
useful for quantum computation. The states that are useful 
for quantum computation should have large state complexity 
and small preparation complexity. 

work of measurement-based quantum computation and 
the circuit model of quantum computation. 

Note that the complex subgroup states presented in 
Sec. |IV E| belongs to the class of stabilizer states. They 
have superpolynomial tree size and can be realized by 
a quantum circuits consist of 0(n^/log n) number of 
gates [H]. Therefore, these states belong to the bottom 
left corner of Fig. But they are not useful for quan¬ 
tum computation since stabilizer circuit can be simulated 
efficiently on a classical computer EH El SHI- 

A. Measurement-based quantum computation 
(MBQC) 

There are several theoretical models of quantum com¬ 
putation, including the circuit model and the MBQC 
model. For the circuit model, the input state can al¬ 
ways be the simple product state. The quantum power 
of the computation lies in the gates applied for coher¬ 
ently manipulating single qubits and entangling differ¬ 
ent qubits [H]. On the contrary, for MBQC, after the 
initial resource state is prepared, we perform projective 
measurements on single qubits and feedforward the re¬ 
sults for choosing the basis of the next round of measure¬ 
ment El m ■ Loosely speaking, all the quantum advan¬ 
tage is contained in the resource state. If this resource 
state is simple, then MBQC will not offer any real speed 
up over classical computation. To make this intuition 
more rigorous, we proved that: 

Theorem 10. If the resource state has TS = poly(n), 
then MBQC can he simulated efficiently with classical 
computation. 

Proof: Consider the resource state in its minimal tree rep¬ 
resentation, one sees that at the lowest layer there are a 
polynomial number of leaves. We will show that it re¬ 
quires polynomial effort to update the tree given a mea¬ 
surement outcome: Assume we measure the zth qubit in 
the basis {I?]), Iry*")} and obtain the the result jr;), then 
for every leaf containing qubit i, say Cq, |a)j -|- C/s |/3)j, we 


update it to (cq {r]\a) -I- cp {r]\/3)) |ry). This requires eval¬ 
uation of the inner products for a polynomial number of 
leaves. The size of the tree after updating can only get 
smaller and thus is still polynomial. So, both the tree rep¬ 
resentation of the state at each step of the computation 
and the update of the state after a measurement can be 
carried out with polynomial effort. It follows that MBQC 
on resource states with polynomial TS can be simulated 
on a classical computer with polynomial overhead. □ 

B. Weaker version of the TreeBQP conjecture 

For the circuit model, rather than checking for each 
algorithm, one would like to have a general proof that 
small tree size does not provide any computational ad¬ 
vantage. In [7], Aaronson raised the question of whether 
TreeBQP = BPP. This remains an open conjecture, here 
we prove a weaker version of it. 

First let us define what TreeBQP is. Bounded-error 
quantum polynomial-time (BQP) is the class of decision 
problems solvable with a quantum Turing machine, with 
at most 1/3 probability of error. TreeBQP is essentially 
BQP with the restriction that at each step of the com¬ 
putation, the state is exponentially close to a state with 
polynomial tree size. In other words, the TS^ of the state 
is polynomial with e = (See Eqn. Since we 

impose more restrictions, clearly TreeBQP C BQP. We 
can also simulate BPP, the classical counterpart of BQP, 
in TreeBQP: One simply implements reversible classical 
computation, applies a Hadamard gate on a single qubit 
and measures in its computational basis to generate ran¬ 
dom bits if needed. Since each classical bit string can 
be represented by a quantum product state, TS is n at 
every steps, so this simulation is in TreeBQP. Thus, we 
have [7]: 

Theorem 11. BPP C TreeBQP C BQP. 

If TreeBQP = BPP, then large tree size is a necessary 
condition for quantum computers to outperform classical 
ones. Unfortunately, we can only prove a weaker version 
of this. For this purpose, we first show a proposition that 
relates tree size and Schmidt rank. 

Note that one can draw a rooted tree in a binary form 
(each gate has only two children) without changing the 
number of leaves (its size). Next, for any gate w we 
denote S{w) as the set of qubits in the state described by 
the subtree with w as the root. Let Y\Z he a, bipartition 
of the qubits into two sets Y and Z. A. ® gate is called 
separating with respect to Y\Z when at least one of its 
children u has the property S{u) C T or S{u) C Z. A 
0 gate is called strictly separating if its children mi,U 2 
satisfy S{ui) C Y and S{u 2 ) C Z. Then, 

Proposition 12. For a bipartition of the qubits into Y 
and Z, if there exists a polynomial sized tree such that all 
the Z) gates are separating with respect to Y\Z, then the 
Schmidt rank of the state with respect to the bipartition 
Y\Z is polynomial. 
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Figure 9. A + gate that joins two (g) gates with the following 
property: One of its children are contained in the set of qubit 
Y and the other contained in Z. The sibling of such 

a + gate must be strictly contained in either V or Z, for (g) 
being separating. Now we can exchange the order of + and 
(g) by distributing \4>Y'2) to \ipYi) and This + gate has 

the same property as before; and this process can be repeated 
upward until it reaches the root, transforming the tree into a 
form similar to the Schmidt decomposition. 


Proof. Identify all the strictly separating (g) gates in the 
binary tree. Since the number of leaves Nr is polynomial 
and the total number of gates in the binary tree is Nq = 
Nl — 1, the number of strictly separating gates, Ns, is 
also polynomial. It is clearer to look at a representative 
example in Fig. Focus on the + gate that joins two 
such (g) gates, \'.pYi) ® Wz) and \^'y^) ® Wz)- Since this 
+ gate contain qubits in both Y and Z, and the (g) gate 
at the top is separating, the qubits under the sibling of 
the + gate must be contained strictly in either Y or Z. 
Without lost of generality, let them be contained in Y 
and denote their state as \4iYf)- We can exchange the + 
gate and the g) gate at the top so that the state becomes 
{\VYf}®\(i)Y.,))®\‘Pz) + {WY^)®\(l>Y.,))®Wz)- Now let us 
relabel \(pYi) ® I'/'Ya) as jlky) and 0 Wy^) as |'I"y.), 
the state can be written as \^Yi)®\vz) and \^'y^)®Wz)- 
The same process can be applied upward until these + 
gates joins at the root. In the final form of the tree, 
one sees that the state has a form similar to the Schmidt 
decomposition: 


N(s 

= ( 45 ) 

Z=1 


criteria, hence do not satisfy Proposition 12 for some bi- 
partitions. For example, the state ( ^—-\ has 

polynomial TS, but there is a bipartition for which the 
Schmidt rank is 


VII. OPEN PROBLEMS 

Even though some properties of tree size have been 
studied, there are still many open problems remained to 
be addressed. Here we list a few of the most interesting. 


A. On the meaning of the tree size 

Is TreeBQP = BPP? If this is true, then the role of 
tree size in quantum computation is clear: Polynomial 
TS means efficient classical simulation, and hence large 
TS is a necessary condition for quantum speed up. 

Is large tree size a resource for any particular task in 
quantum information? 

Given an explicit family of states, can Raz’s theorem be 
modified to provide an exponential lower bound, instead 
of 

Is there an algorithm (other than exhaustive search) 
to find the optimal tree given a quantum state? 

How to prove and verify the superpolynomial TS of 
mixed states? 

Has any (family of) states with superpolynomial tree 
size been produced in experiments? 

As mentioned in Sec. 3, the ground state of a ID 
gapped Hamiltonian, which can be described by an MPS, 
has polynomial TS. Is there a physically reasonable ID 
two-local Hamiltonian whose ground state has superpoly¬ 
nomial TS? A ground state at phase transition no longer 
obeys the area law because of high entanglement. At this 
point the state is not an MPS so one can expect that its 
TS is large. 


B. Technical open problems 


where |'I'y)j contain qubits in Y and z)i qubits in Z. 
Ns, the number of terms in this Schmidt-like decompo¬ 
sition upper bounds the true Schmidt rank, hence the 
Schmidt rank is polynomial. □ 

Now suppose that at every step of the quantum compu¬ 
tation, PropositioifT^ is satisfied for all bipartitions, then 
the Schmidt rank is polynomial for all bipartitions. It fol¬ 
lows from a theorem by Vidal [3] that the computation 
can be efficiently simulated with classical computers. 

There are states that do not satisfy the condition of 
Proposition one example is the optimal tree of the 
most complex four qubit states (see Eqn. ©)• There are 
also states with polynomial TS that do not satisfy Vidal’s 


Below are a few more technical open problems: 
e. Tree size 2" for n qubit states: One observation 
we made for the tree size of a few qubits is that the most 
complex state has TS = 2" for n = 2,3, 4. Is this a pure 
coincidence, or is this generally true for any n? 

/. Stable tree size: Eor the three and four qubit case, 
the most complex state are unstable. Infinitesimal per¬ 
turbation in suitable directions in the Hilbert space could 
reduce its tree size to the second most complex class, 
hence the maximal value of the stable tree size TS^ is 
different from the maximal tree size. Is the maximal sta¬ 
ble tree size always smaller than the maximal tree size? 
Is it always equal to the second largest tree size? 
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g. Exponential tree size: The TS of the permanent 
state and the 2D cluster state are shown to be superpoly¬ 
nomial. We conjecture that they in fact have exponen¬ 
tial tree size, TS = and TS = with some 

0 < e < 1/2 respectively. There are strong evidences to 
believe these two states have exponential TS: for the per¬ 
manent state it is known that computing the permanent 
of a matrix is ^P-hard |36j . and a subexponential tree 
size for a permanent state would imply a subexponential 
formula to compute permanent, contradicting the expo¬ 
nential time hypothesis (or a variant of it, ^ETH) [? ]; 
for the 2D cluster state, if we assume the contrary that 
it has subexponential tree size, by Theorem we can 
simulate the polynomial-time quantum factoring proto¬ 
col with some subexponential effort. This contradicts 
with the belief that quantum computing offers e^onen- 
tial speed up compared with classical computing Since 
2D cluster states are 2D PEPS with a bond dimension 
X = 2 [? ], so TS(|2D)) = 20(v^) and thus e is upper 
bounded by 1/2. 

VIII. CONCLUSION 

In this paper, we revisited the complexity measure, 
called tree size, for pure n-qubit states. For few qubits, 
the state with the largest tree size is identified for n = 
2,3,4 qubits. For 4 qubits, the most compact states ad¬ 
mits an optimal tree that is not recursive. The general¬ 
ization of tree size to incorporate small fluctuation and 
to mixed state is also discussed. 

Raz’s theorem on the superpolynomial lower bound of 
multilinear formula size can be utilized to show that some 
multiqubit states have superpolynomial tree size. Exam¬ 
ples of such complex states, the Immanant states and 
the subgroup states, are described. Moreover, the con¬ 
jecture that the 2D cluster state has superpolynomial TS 
is proved. We also show how to verify the superpolyno¬ 
mial TS of stabilizer states, such as the complex subgroup 
states and the 2D cluster state, with polynomial effort by 
measuring a stabilizer witness. 

The relation between tree size and quantum computa¬ 
tion is discussed. In measurement-based quantum com¬ 


putation, if the initial resource state has polynomial tree 
size, then the computation can be simulated on a classi¬ 
cal computer with polynomial overhead. For the circuit 
model of quantum computation, we show that most of 
the states arising in Deutsch-Jozsa algorithm have large 
tree size. A similar result for Shor’s algorithm is also re¬ 
viewed, although a number-theoretic conjecture need to 
be made in this case. Finally, we present a proof for a 
weaker version of the TreeBQP = BPP conjecture, which 
says that if the tree size of the quantum state is polyno¬ 
mial through out the computation and obeys some extra 
conditions, then the computation can be simulated effi¬ 
ciently. 

In conclusion, tree size of quantum states as a com¬ 
plexity measure possesses some desirable properties. The 
most important one is that it is possible to derive non¬ 
trivial lower bound on tree size. We have seen some signs 
on the complex relation between tree size and the use¬ 
fulness of a state for quantum computation speed up, 
but the picture is still unclear. By further investigat¬ 
ing tree size and other complexity measures, we hope 
to identify state complexity as the resource for quan¬ 
tum computation. With that understanding one can rule 
out states which do not provide quantum advantage, and 
concentrate on producing and characterizing states with 
high complexity, and possibly identify new quantum al¬ 
gorithms based on complex states. 
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